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Abstract. This work develops the asymptotic properties (weak consistency 
and Gaussianity), in the high-frequency limit, of approximate maximum likeli- 
hood estimators for the spectral parameters of Gaussian and isotropic spherical 
random fields. The procedure we used exploits the so-called mexican needlet 
construction by Geller and Mayeli in |21l . Furthermore, we propose a plug-in 
procedure to optimize the precision of the estimators in terms of asymptotic 
variance. 



1. Introduction 

The aim of this paper is to investigate the asymptotic behaviour of a Whittle- 
like approximate maximum likelihood estimates of the spectral parameters (e.g., 
the spectral index) of isotropic Gaussian random fields defined on the unit sphere 
S 2 . We employ a procedure based on the so-called mexican needlet construction by 
Geller and Mayeli in |21j . Furthermore, we develop a plug-in procedure aimed to 
merge and to optimize these results with the achievements pursued in [12] , [13) , see 
also [14] . where the asymptotic behaviour of Whittle- like estimates were studied 
respectively in the harmonic and standard needlet analysis frameworks. 

Under the hypothesis of Gaussianity, fixing smoothness conditions on the be- 
haviour of the angular power spectrum, we pursue weak consistency and central 
limit theorem allowing for feasible inference. From the technical point of view, the 
asymptotic framework we use here is rather different from the usual, being based 
on observations collected at higher and higher frequencies on a fixed-domain (i. e. 
the unit sphere) . In this sense, this work can be related to the area of fixed-domain 
asymptotics (see for instance [2], [34]); on the other hand, as for [12] and |13j . some 
of the techniques used here are close to those adopted by [46j to analyze the as- 
ymptotic behaviour of the semiparametric estimates of the long memory parameter 
for the covariance of stationary processes. In terms of the angular power spectrum, 
we shall also focus on semiparametric models where only the high-frequency/small- 
scale behaviour of the random field is constrained. In particular, we consider both 
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full-band and narrow-band estimates, where the latter allow unbiased estimation 
under more general assumption, by paying the price of a slower rate of convergence 
if compared to the former. 

This investigation, as many others regarding statistical inference on spherical 
random helds, is strongly motivated by practical applications, especially in cosmol- 
ogy and astrophysics (see for instance [3TJ and the references therein). For instance, 
as described in [9] and [8], satellite missions such as WMAP and Planck are now 
providing huge datasets on Cosmic Microwave Background (CMB) radiation, usu- 
ally assumed to be a realization of an isotropic, Gaussian spherical random held: 
the issues concerning parameter estimation have been considered by many applied 
papers (see [24] . |31j for a review), but in our knowledge, until now, rigorous as- 
ymptotic results are still missing in literature. We however refer also to [4], [15j . 
|19j . [43j . [44j . [36j for further theoretical and applied results on angular power 
spectrum estimation in nonparametric settings, and to [25], |27j . [26 , [28 , 32J, 
[29j and [37] for further results on statistical inference for spherical random helds 
or wavelets applied to CMB radiation. 

Another result we work out in this paper concerns the formulation of a plug-in 
procedure which combines the application of the asymptotic results here attained 
with those described in 1121 and [13] , where the authors proved that weak consis- 
tency and central limit theorem can be achieved respectively by standard Fourier 
and standard spherical needlet analysis. In |12) . the authors themselves have put 
in evidence that, if the asymptotic achievements are better with respect to those 
obtain in needlet framework in terms of precision of the estimates (e.g. their asymp- 
totic variance is smaller), in many practical circumstances the implementation of 
spherical harmonics estimates may present some difficulties, due to their lack of lo- 
calization in real space. The presence of unobserved regions on the sphere (common 
situation in the case of Cosmological applications), can indeed make their imple- 
mentation infeasible, and spherical harmonics exclude the possibility of separate 
estimation on different hemispheres, as considered for instance by [5], [45] . In view 
of these issues, in |13j . the authors investigated the Whittle- like procedures to a 
spherical wavelet framework, in order to exploit the double-localization properties 
(in real and harmonic space) of such constructions, at the cost of a smaller precision 
in term of convergence in law of the estimates. They focussed their attention on 
spherical needlets, second-generation wavelets on the sphere, introduced in 2006 by 
[40) and [41] . and very extensively exploited both in the statistical literature and 
for astrophysical applications in the last few years: for instance, their stochastic 
properties are developed in [4], [5], [6] [29] . [30] and [39] . More recently, needlets 
have been generalized in different ways: we cite spin needlets (see |17j ). and mixed 
needlets (cfr. |18j ). which represent the natural generalization to the case of spin 
fiber bundles, again developed in view of Cosmological applications such as weak 
gravitational lensing and the polarization of the Cosmic Microwave Background 
(CMB) radiation (see for instance 0], [S], p], OS* PU* [IS], [38], SU, [45] . 
[47]). On the other hand, needlets have been generalized to an unbounded support 
in the frequency domain by [20] . [21] and [22] . the so-called Mexican needlets. 
In this case, as we will describe in details below, even if the support in frequency 
domain is unbounded, the form of the weight function, depending on the scale pa- 
rameter p, is such that for each wavelet there is a small numbers of frequencies 
which give a contribution substantially far from zero , while in the real domain the 
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same weight function allows a closer localization than the one related to standard 
spherical needlets. In particular this double localization depends on the value of p 
or, better, on its distance from the spectral index, allowing these estimates to be 
more efficient than the ones obtained with standard needlets. Our idea, therefore, 
is to build a plug-in procedure on two steps, the first step being to estimate ap- 
proximately the value of the spectral index by standard needlets and the second 
step providing a estimation with mexican needlets, whereas the value of the scale 
parameter p will allow a more efficent estimator. 

The plan of the paper is as follows: in Section^ we will recall some background 
material on mexican needlet analysis for spherical isotropic random fields; in Section 
[3] we will introduce and describe the Whittle-like minimum contrast estimators, 
while in Section [4] we shall establish the asymptotic results on these estimators. 
In Section [5] we present results on narrow band estimates, while in Section [6] we 
will describe the plug-in procedure mentioned above. Finally, the appendix collects 
some analytical and statistical auxiliary results. 

2. Random fields and mexican needlets 

In this Section we will introduce the mexican needlet framework (for more 
details, cfr. |21j ) and its application to the study of the isotropic, Gaussian ran- 
dom fields on the sphere. First of all, consider the set of spherical harmonics 
{Yim '■ I > 0, m = —I, /}. As well-known, it represents an orthonormal basis for 
the class of square- integrable functions on the unit sphere space L 2 (§ 2 ): the spher- 
ical harmonics are defined as the eigenfunctions of the spherical Laplacian A 52 cor- 
responding to eigenvalues — 1(1 + 1) (see, for more details and analytic expressions, 
PQ [49j . [50] . [37] and, for extensions, [33], [35] ). The mexican needlets are defined 
in [21] as 

(2.1) V ifc;p (*) := fp E F '™ (*) F <™ fefc) > 

l>l ^ ' m=-l 

where 

(2.2) f p (x) = x 2p exp (-x 2 ) . 

Observe that is a set of cubature points on the sphere, indexed by resolution 

level index j and the cardinality of the point over the fixed resolution level k, while 
Xjk > corresponds to the weight associated to any £ - fc . The scalar Nj denotes 
the number of cubature points for a given level j (cfr. [40] . [41] , see also e.g. [21] 
and [37] ) , chosen to satisfy the following 

(2.3) \ jk « B~ 2 i , Nj « B 2 i , 

where by a sa b, we mean that there exists cx,C2 > such that c%a < b < C2C1. 
Below, we shall assume for notational simplicity, as in [13] . that there exists a 
positive constant cb such that Nj = cbB 2 ^ for all resolution levels j. In practice, 
cubature points and weights can be identified with those evaluated by common 
packages such as HealPix (see for instance [4], [10] . [23 ). 

Considering L t ((x,y)) = J2 m =-i Yim ( x ) Yim (v) as a projection operator, the 
definition (|2.1[) corresponds to a weighted convolution with a weight function (|2.2p : 
mexican needlets can be considered as an extension of the spherical standard 
needlets, proposed in [40] . [41j . see also [6], [12] . [37] . The main difference between 
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these two kinds of wavelets concerns their dependence on frequencies. In fact while 
standard needlets have a compact frequency support (see again [40j . [41]). each 
mexican needlet is defined on the whole frequency range. In [211 . mexican needlets 
are proved to form a nearly tight frame, differently from the standard needlets 
which describe a tight frame and, as consequence, are characterized by an exact 
reconstruction formula (see again |40j ) . 

Consider now a zero-mean, isotropic Gaussian random fields T : § 2 x SI — >• M: 
it is a well known fact that for every g £ SO (3) and x £ S 2 , a field T (•) is isotropic 
if and only if 

T{x)=T (gx) , 

where the equality holds in the sense of processes (see [36], [37] ), and that (see e.g. 
|37j ) the following spectral representation holds: 

T ( x ) = a hnY[ m (x) , a lm = T (x) Y lm (x) dx . 

l>0 m— — l 

Note that this equality holds in both L 2 (§ 2 x fi, dx ® Pj and L 2 (P) senses for every 
fixed x E Ei 2 . For an isotropic Gaussian field, the spherical harmonics coefficients 
aim are Gaussian complex random variables such that 

E(a; m ) = , E(aj m a/ imi ) = S^S^Q . 

The angular power spectrum {Ci , I — 1, 2, 3, ...} fully characterizes the dependence 
structure under Gaussianity. Properties of the spherical harmonics coefficients un- 
der Gaussianity and isotropy are discussed for instance by [3], [37] ; here we recall 
that 

i 

^2 I '™' 2 ~ Cl X *2Z+1 • 
m— — / 

Hence, given a realization of the random field, an estimator of the angular power 
spectrum can be defined as: 

1 1 

m——l 

the empirical angular power spectrum. It is immediately observed that 

(2.4) E( a <) = 2lTI ± a=C,,^(§)=^^0for^ + oc. 

in— — I \ / 

As in [13] , we introduce the following regularity condition on the angular power 
spectrum: 

Condition 1 (Regularity). The random field T(x) is Gaussian and isotropic 
with angular power spectrum Ci so that for all B > 1, there exist ao > 2, cq > 
such that: 

(2.5) Ci = r ao G (I) > 0, for allien , 

where Cq 1 < G (I) < cq for all I € N , and for every r G N. there exists c r > such 
that: 

d 



d U r G ^ 



<c r u r , e (0, +oo) 
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This assumption is fulfilled by popular physical models, for instance in a CMB 
framework the Sachs- Wolfe power spectrum, which is the leading model for fluctua- 
tions of the primordial gravitational potential, takes the form (I2.5[) . see for instance 

m 

First of all, we stress that Condition Q] implies the following Condition [2l given 

in 



Condition 2. Condition [JJ holds and, moreover, there exist ctQ > 2 and a 
sequence of functions {gj (•)},-=! 2 such that: 

(2.6) d = r a ° 9j (^] > 0, for all B^ 1 < I < B* +1 , j = 1,2... 



B3 

-1 



where c < gj < cq for all j G N , and for every r — 0, Q, Q £ N, there exists 
c r > such that: 

d r 



sup sup 

j Bi- 1 <u<Bi+ 1 



du 



~9j (u) 



< Of. 



As an example, consider 

Cl ~ l m; 

where P (I) — Yli=i c vA % an d Q (0 = Si=i c qM l are positive polynomials of degree 
p and q respectively, so that ao = a — p + q > 2., so that 

^ _J_ C P.P~1 _l _J_ c p,p-2 1 _|_ 



1 



-1 1 I C q,q-2 1 



-q,q 1 -1 " T 1 J 77 "T •■■ 

2 



7— a+p— q °P,P ''''' ' i> 



1 Cp^p-iB 3 j 1 Cp.p-2 / B J 



Cq.q 1 j 1 Cg.g-i 1 Bi , 1 c g.<;-2 1 / Bi V 



- rao n^ 

Condition [1] will be necessary to prove needlet coefficients (|2.7[) to be asymp- 
totically uncorrelated (see |30j . [39]); as we shall show, Condition [1] is sufficient 
to establish consistency for estimator we are going to define but we will consider 
two further nested assumptions, [3] (which implies and is implied by [T|), to obtain 
asymptotic Gaussianity, and 2] (which implies [3]) to provide a centered limiting 
distribution, see also 



Condition 3. Condition^ holds and moreover 

G(l) = G (1 + nl- 1 + O (I- 2 )) . 
Condition 4. Condition^ holds and moreover 

G(/) = Go(l + o(r 1 )) . 
For any given j, k, p, we define the needlet coefficients as: 

Pjk;p : = / T 0) ^jk;p 0) dx 



s 

/ 7 \ * 

2>1 v 7 m=-l 



(2.7) = v^E/pI^J E^^fe*) . 
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SO that 

E (P jk . p ) = v^^/p ( ) E E ^ Y ^ fefc) = 



I - 1 



2>1 x ' m=-i 

Under Condition [51 the following result is given in |30j and 

Lemma 2.1. J/ < 4p + 2 — ojo < Q, then under Condition^ there exists a 
constant Cq > 0, such that 

Corr (p jk . p ,p jlk ,. p ) < 



[i + s((i+i')/2-io gB [a+i')/2])d (^ k ,^, k ,)] 



(ip+2-ao) ■ 



Assume now that from the observations over the random field, we are able to 
build the following set of quantities 

/ I \ ~ N] 

:= E # ( bJ ) (2/ + 1} Ci ~ E for each J ' G t Jo ' Ji l ' 

i>i ^ ' fe=i 

where the last approximation is motivated by the nearly tight frame property, as 
in [39] . 

The next result describes the asymptotic behaviour of the variance-covariance 
matrix of Aj- P in terms of j. 

Lemma 2.2. If Condition^ holds with < 4p + 2 — a a < Q, fixed Aj £ Z, we 
have 

1 /— \ 2G 2 

^ B 2(l- aoh Var { A r,p) = 4 4 P+ (l-a ) r (4P + 1 - «o) ; 

/™ ^(T^)7 Cou (AiiP. A ^a j; p) = 2Gg 3 ^M y L (4p + 1 - a ) , 

(2.8) r p (Aj) := B A J'(i-«o) C osh (Aj log B)- (4p " ao+1) . 

Proof. Simple calculations lead to: 

Var (A j;p ) = Var \J^P%, P ) = ]T f A p (-^) (21 + if Var (c x 



Jfe=i / i>\ 



B3 

1>1 x 



while, for Aj £ Z, 



Co«(A JiP> A i+Aia ,) = C i;[ £>2* liP , £ ? 



j+Ajk 2 :p 



,fel=l fe 2 =l 
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= c ™ E fl (w) + 1) c h , E {^kj ) (2fa + 1) a 

\ii>i v 7 ; 2 >i v 



2 E^wk 2 (^)^ + 1 ^ 2 



Under Condition [TJ by applying Lemma IA.21 in view of the equation (IA.1[) with 
a = 4 and n = 1 — 2ao, we have: 



= 2G 



£2(l-a )j 



I -^ T r(4p + i- ao ) + 0j (s 2 ( 1 -^) , 



r 4 4p+(l-c 

while, for the equation (IA.2|) with a\ = a% = 2, n = 1 — 2«o and t p (Aj) — 
t p,2,2 (Aj), we obtain: 

COV (Aj ;p ,Aj + Aj;p) 

- <E/, 2 (^)4 2 (^)^ 2Q0 +^D 

= 2^ 4 y +( ( 1 A £ ) i3 2 ( 1 -°^r (4p + 1 - a ) + o (s^"-) 



as claimed. □ 

3. Mexican Needlet Whittle-like approximation to likelihood function 

In this Section, our aim is to define a mexican needlet Whittle-like approxima- 
tion to the log-likelihood function of isotropic and Gaussian random fields on the 
unit sphere under Condition [T] and to develop the corresponding estimators. We 
will follow a strategy analogue to the one used by [13] . (see also [12] and [46] ). We 
let 

where f3j k . p is defined as in (|2.7[) . Again, under the hypothesis of isotropy and 
Gaussianity for T, we have 

where 

r= [Cov(p jk;p ,p jk , p )] Kk , = ±- |e/ 2 (2/ + i)aj i Nj , 

in view of (|2.7[) and Lemma [2.11 (see also |12) . [13]). The likelihood function is 
then defined as 

£ (": •/,/.) = (2tt)-^ (detL)- 1/2 exp|-^j; p L- 1 t J;?) | • 

Let 
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K 



i>i 



Under Condition [TJ we have: 

C [a, G- ~$ 3 , p ) = (2tt)-"' {GKf (a))"^ /2 exp ■ 

and the corresponding approximate log-likelihood is 

-21og£(a,G;^ > )=X:|^- y -log 

up to an additive constant. 

By summing with respect to j, we obtain. 

-l 



1 12k Pjk;p 

2 GKf (a) 



H 2 



GKf 1 (a) 



Jl 



Jl 



Rf ojL (a, G) :=\J2 N i £ ~ 2 lo g £ ( a < G > P. 

\j = Ja J j = Jo 



where the choice for Jq, Jl will be discussed later. Hence we define (cfr. [12] and 



arg mm Rf o JL (a,G) , 



(a,G)ee 

where O = [2, +oo) x (0, +oo) . Computing the derivative of Rj* j with respect to 
G and setting it equal to zero, we have 



L>j=Jo iV J j=Jo 



dG 
whence 

Since 

Jl 



J2k P 2 jk; P Nj 



G 2 Kf{a) ' G 



J2 Jl t Ni f-f Kf (a) ~ J2 Jl j N< f-i Kf (a) ' 

^3 = Ja J J = Jo J ^J = Jo ■> 3 = Ja ■> 



Jl A- 



t- 



d 2 

qq5 R To,Jl ( a ' G ) 
2A- 



*-'3=Ja 3 3=Jo \ ■> 



G=G¥ , (a) 
■'a • j l v ' 



(G%,j L («) 



> 



and j (a,G) — >• +oo as G -> or oo, the second derivative test yields that 
Rj^ Jl (a, G) has a unique minimum over the domain on Gj^ Jl (a) . Therefore, we 



Jo,Jl 

can define 

M 



i? 



1 + logG^ 



tff (a) 
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Remark 3.1. In this formula it is necessary to fix explicitly the values of L, Jq 
and Jl- Let us fix L as the highest multipolc level available from the dataset. Given 
L, as stressed above, differently from the standard needlet case (see for instance 
40 . [41 ), in the mexican needlet case the weight function does not have a compact 
support. Therefore, for computational reasons, we must find a criterion to choose 
the suitable extrema for the sums over j involved. Considering (see again [21] ) the 
behaviour of f p (•), we can fix thresholds Sb,i (L), Eb,2 (L), such that: 

Jo = max | j el: f p (jjj+r) > £b,i (L) f p (Jj- 
J L = min j j e Z : f p ( — ) < e B ,2 (L) f p 



If, for instance, we choose, 



SB.i (L) = ^_ cxp , cb,2 (L) = ^_ exp (i? 2 (i? 2 - 1)) 

we find B J ° = B, B Jl = L/B, similarly to the classical needlet case as described 
in [15] . 



4. Asymptotic Properties 

In this Section, we prove weak consistency for the estimators aj' Jl and G^f a Jl , 
and for the former also asymptotic Gaussianity. We begin with some definitions: 
let 

1 ^ „ G Kf (ao) 



Gj ;J L («) ~ e j l E K M {a) ■ 

^3 = Jo J ] = Jo 3 v ' 

Computing the first and second order derivatives of G I j I o . Jl (a), indexed by n, we 
obtain 



G jijL,n (a) ■ = -^ G Jo;Jl («) 



Go ^ M Kf*(ao) TT 



■>j=Jo 3 j = Jo 

where (see Proposition [T]) in the Appendix, we have 



(4.1) r /:! (ni : ( - T -^f— ) = (jlogB + ^— + o l{ \) 

2 



Uy,i (a) = 2 



Jiff (a) V K™(a) 
Kf(a)) Kf(a) 

2 



(4.2) = f log 2 B + 2j log^flM + 2 - + 0i (1) 

Furthermore, we fix 

t/ J - ; o(a) = l,G^ (a) = G^ Jt (a) , 
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(since now, we will use either G 1 f Q .j L (a) or G 1 j q . Jl (a)). Recalling that Nj = 
C B B 2 K Thus by JO), we have for's = 0, 1, 2, 

^j=Jo a j=Jo 

The next auxiliary result is as follows: 

Lemma 4.1. Assume Conditional] holds with < 4p + 2 — ao < Q. M^e /icwe 



lim E 

Jl,— too 



{GjU (ao)) G , 



w/iere 

2 r(4p+l-a ) 



°o : = a o (p,"o) 



' 2 4p-ao p2 (2p - ^ + l) ' 

and To is as defined in Lemma \B.l\ 
Proof. We have 

■lL E ( A : 

E ' 



J2 Jl , Ni ^ T Kf (a) 

^3 = Jo ■> J = Jo 

Go A E ; ft (jr) (2/ + 1) l~ a ° (1 + (r *)) 
Y Jl t ^ Kf (a) 

^3 = Jo 3 3 = Jo 3 v ' 

= Go + o Jl (1) . 
On the other hand, we prove that 

(4-3) Cov ( r A » Ai » ) = c\olB^B^r B (Aj) . 

\G Kf (a) G K^_ Aj (a) I 



We can indeed observe from Theorem 12 . 2 1 that 



Cov 



G Kf(ayG K™ Aj (a) 



2 T2 („\ U-2 a , C ° V { A 3^3+A3) 



B aAj 2r(4p+l-a ) 
Il (a) ^v+i- 
4pH 
(27 



-TB 



(Aj) B 2 ( 1 + a - a o)3 



T- I , > I 44p+l-Ct(: 

(4.4) = 24 r (4p + 1 - Og) ^aAj^d+a-ao),- 

1 ; 2 4 p- 2q «+« r 2 (2p - % + 1) v ; 
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Hence 



Var 



Go 



1 ( Jl A- Jl ~ j 



A 



j+Aj;p 



1 V- si ( Sfci ^jfcljP Efc 2 Pj+Aj,k;P \ 

i ^ W («) ' Go*** («) J 
1 1 r(4p+l-^o) ^ B 2(l + a- Qo)j A , 

(ES, ^) 2 42P " a0+fr2(2p " f + 1) ^ ^ 



Following Lemmas IB . 1 1 and [2| and computing in a — oto, we have 



Var 



Go 



2(1+ to) r(4p+l-a ) 



5 2 - 1 



2 




_B 2 



al(l + r )B-^+o(B-^) 



Lemma 4.2. Under Condition]]^ we have for s — 0, 1, 2: 



□ 



sup 



G 3U;» 



Proof. Under Condition [TJ we can readily obtain that 



7Zl~ L 7 T G a Kf( ao ) 77 



/3 



,G Kf (a ) 
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SO that 



( P%;p 
\G Kf( ao ) 



( 



< 



(Jl + Jo + 1) 



V 



><5 e 



X (J L + Jo + 1) 

In view of (|A.4[) and (IA.5I) , we obtain 



>S E 



E 



E^ L j B^ 1+a ^ a °^ j s 

g(2+a-a ) _ I J L J^JL(l+a-a ) _ J ^Jo(l+Q-ao)-l 



so that 



sup 

j 



{J L + Jo + 1) 



- l) J L B jL< - 2+a - a °) - J B J o(2+a-Q )-l 



^E 



a 



< +0O 



On the other hand, we have by Chebyshev's inequality and Lemma l4.il that 



< 



1 



Goi^f (a ) 



> <S e ( J L + Jo + 1)' 



£ 2 ( Jl + Jo + I) 2 



Var 



^E 

/a7 v 



jfc;p 



Go AT (a„) 



- 1 



whence 



sup 

\3 = Jo,-Jl 



{J L + Jo + 1) 



$ jk;p 



< ( Jl + Jo + 1) sup 1 

3=J ,---Jl 



G Kf (a ) 



^E 



B 2 



3 k 



Go A" (a ) 



> <5 e ( J L + Jo + if I 

> <J e (Jl + Jo + I) 2 
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< 



(J L + Jo + 1) 



□ 



Let us focus now our attention on consistency, following a technique developed 
in [7] and used in [46j for long memory processes (see also |12j and |13j ). 

Theorem 4.3. Assume Condition^ holds with < 4p + 2 — ao < Q, we have, 
as Ji ^ oo . 

a J ,J L p q o > 

Proof. Following |46j (see also [13] for the standard needlet case), we let 
AR% jL (a,a ) = R^ Jl (a) - B% tJt (a ) 

g jLj l («) 8 Go 
^ G^>) 1 ^ Kfja) 

= ^(a)-T^(a) , 

where 

Gf j (a) Gf , (a ) 

(4.5) r« JlM = ,o gj ^-l„ g -^, 

Gf r (a) 1 A iff (a) 

(4.6) = log ^ii + £ Wjlog _^. 

For any e > 0, we have 

P ( S 1.7, "«o >e) = Pf, mm AijJ^ (a, a„) < ' 

= pf min T% tJli (a) + U% tJL (a)<0 

\\a — ao|>£ 

Hence, by combining Lemmas 14.41 and 14.51 we obtain 

I™ U"j L (a,a o )>0, 

as claimed. □ 
Lemma 4.4. Let Uj 1 oJl (a, ao) be defined as in |^.6'[ j. For all e < ccq — a < 2 

lim Uf j (a, a ) 

Jl—^+oo 



lim log — -j > , 7 > Nh log ; ; 
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B 2 — 1 B 2 
log Bi2+a - ao ) Z i + lo Z B B^T a ~ Q ° > * e > ° ' 



&/ «o — a = 2 we have 



lim 



log Jl 



Uj ,j l (a, a ) = 1 



and if ao — a > 2 we have 
lim 



1 tt I \ a ° ~ 01 1 
-UJoJl («.ao) = — ^ 1 



log B J t 

Proof. Consider first the case e < ao — a < 2. For the sake of simplicity, we 
fix Jo = — Jl. We have that 



■'3=—Jl ■? ]=-Jl 

= (a- a ) \ogB ^J L 
On the other hand, we have 



— - ]T N i ( lo S B^- a ^I p {B, a - ao ) + o(j)) 
33;) + log (ip (-B, a - a )) + OJi (f ) 



log — V N- 



W £ B 2j . ^ J iVf ( a ) 

^>3=—Jl j=-Jl 
Jl 

ao)3 +oj L {l) 



I p {B,a-a ) y B2jB(a 



= l0 S B 2+(a-ao)_l g(a " a0)(Jl, + 1) + l0 § ^ ^' ° - a °)) + °^ W 

_B 2 — f 

= log B 2 +(a - ao ) _ I + ( a " a o) ( j l + 1) logS + log (I p (S, a - a )) 
As shown in [13j, we have that the function 



l(x) 



B 2 -l ( B 2 log B 

B 2 + x -1 +X ( B 2 - 1 



has a unique minimum at x = 0. Therefore, for any \a — ao\ > e > 0, there exists 
a constant <5 e > 0, such that 



If a — ao < —2, we have 



Uj ,j l (a,a ) > S s . 

1 , 



1 



log B 2J 

ao — a 



log 



Jl 



B j(2+a-a ) 



n, (a— ao) , T , Kf 1 (ao) 
log S 2Jl - ^- 2^ AT,- log ' 



i=Jo ■? j=Jo 



Kf{a) 



- 1 



o.j L (1) 
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Finally, we have for a — ao = — 2 



lim — 

Ji->oo log J L 



Uj ,j l (a,ao) 



lim — 

Jt-i-oo log Jl 



{-logB 2J - + log J L +0 Jl (1) + log B 2Jl 0j l (1)} = 1 



□ 



Lemma 4.5. As Jl — > +oo, we have 

su p\ t j^,j l ( a > a o)\ = o p (1) . 

a 

Proof. Because 



'j=Jo 3 Kf'(a) 



it follows from Lemma 14.11 that 



while 



Indeed, we have 



Var 



(Gf j (a) \ . N 



(a) 



_ (j^+a-ao) _ !) 2 /?)(j3 , a - Qo) 2 r (4p+l-a ) 2J , f aJj; , 

(ij2(l+a-ao) _ iJao-a42p-a + f p2 fop — % + l) ^ 7 

By Chebyshev's inequality we have 

Gf ,j L («) 



^SU («) 



and from Slutsky's Lemma 



loe 



On the other hand, by Lemma 



sup 



as we claimed. 



G%,J L to 



l->,0, 



- 1 ->« o . 



<3U («) 



- 1 



>p , 



□ 



Our purpose now is to study an asymptotic convergence of estimator ctjg j L 
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Theorem 4.6. Let < Ap — ao < Q. Assume Condition^ holds with . Hence 
we have 



oo . 



(4.7) B Jl (a A / - a ) = O p (1) , as J L 

Under Condition^ we have 

(4 - 8) y^j a ) - Ipo(ao) {B + l) > 

Under Condition ^ we have 

(4.9) B Jl (af -a ) -+ d N (0,$) 

where 

(B 2 — l) 3 

«o : = {P, B ,<*o) = o-l (1 +r) 



2 2 r(4p+l-a 

= ct (p,a ) 



2 4P- Q0 (r(2p+l-a /2)) 2 ' 

— ({B 2 + 1) (t + T 2 + ?0? 2 ) + 2?! - T?) 



wii/i To Ti T2 as defined in Lemma \B.l[ 

Proof. Again we shall focus on the Taylor expansion 
d 



= ^RjU («) lo=a5J, j£ = ^o,/, («) la=a„ + Qf ,J L («) |c=* - a) 



da 
da 2 



(«) = («) ; 



where a € [ao — ao + 5j L ] , and 5j — > p as Jl — > oo by Lemma [4.31 The 
equation above then leads to 

The proof is readily completed by combining Lemma 14.71 and 14.81 □ 



Lemma 4.7. Assume Condition^ holds with < 4p + 2 — ao < M , we have 
B J ^S M (a..) ) 4,o(<*o + l) log^ . 

?/ Condition [^] ZioWs we ftave 

B JL Sf 0;jL (a ) iV (o,^ (1+7) 
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Proof. Note, first of all, that, as in [13] . the proof of (|4.8|) is totally equivalent 
to the case of (|4.7|) . First of all, we can rewrite Sj 1 Jl (a) as follows. 



^j=Jo •? j = Jo k 3 

Gl./,,i («) 1 A («) 



Jl jy. A—' 

j=J J j = Jo fe 

Jq,Ji.,i y' x > 1 y^ -^v 

We can easily see that 

^ M = e^K)^'^ (ao) ' 

where 

-qM 1 A K% (ap) v / /3* fc;p gg^ (ap) ' 

^ W ^^2.^(00)2.^0^(00) Go 

and from Lemma 14.21 we have 



GJU.ifo) " 

Under Condition [3] we have 

lim S J -E(^„7 L (ao)) 

e£ Jo ^ v *f ^)\ ^ e£ Jo ^ i G °^f («,,) J 

,. 7p,o(a + l) logB 

= hm - K -L- — +0J (l) ; 

while under Condition [4] we find immediately 

In order to compute the variance of j (oto), we split it into 3 terms (see again 



1 



Var (sf OtJL (a ))=A + B + C , 

where 



18 



CLAUDIO DURASTANTI AND XIAOHONG LAN 



By fixing j = ji, Aj = j 2 - ji, we have: 

1 V yf /^J^o^jWM^ 



Efcj Pjki ;p Efc 2 Pj+Ajk 2 ;p \ 



xCov 



G Kf (aoV G Kf +A] {a ) 



(B 2 -l)log 2 B 2 / t2 ^ 2 



+ B2 + 1 2 (l + r 2 )) +o(B- 2J -) 

(b 2 — i) y V 



by applying Lemmas IB. II and IA.2I On the other hand, from Lemma 14. 11 we obtain 

b = ^EE f ^ (^h (ao) ) 



31 J2 



Ji \ / Jl \ Jl Jl-3 

j I E log B»B 2 » E l°Z&*B 2 » E E B a ^r B (Aj) 



E/=J n ) \ii= J o J \h = Jo J h=Jo Aj=-J -j 

\ 2 



£? 2 — l) log 2 B , / / 1 \ 2 „ ' 

a 2 B- 2Ji Ji-—— (l + r ) 



B 2 " \\ B 2 -l 

Finally, we have that 



E 5 2j2 logB 32 ] +o(B 2Ji 



x (A-BrrT))+°(B 2 ''': 
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Summing all these terms, we obtain 

, B 2 log 2 B 



A + B + C = o% (1 + r) B- 2J - + o {B 2J - ] 



We can use the Lemma \B. II to observe that 

a 2 (1+7) 



WU(*°)) = /°; 1+r \, (^+^(i))- 



Far 



Hence we have 

a 2 (1 + r) B 2 log 2 B 



lim B 2jL Var 



( S 1,Jl (»o) 



j l ^oo V Jo < Ji v U V (B 2 - 1 



To prove (|4.9[) . it remains to study the behaviour the fourth order cumulants, 
observing that this statistics belong to the second order Wiener chaos with respect 
to a Gaussian white noise random measure (see |42) . |13j ). Let 

B J -Sj L {^) = ^j I Y.^ +B ^ > 

3 

where 

f V P 2 

(4.10) Aj = B 2 nogB J 



J2k Pjk 



(4.11) B j = B 2 nogB j 



NjGoKj (o ) 
^(ao) J 



Go 



In the Appendix, Lemma IC . 1 1 proves that 
1 

^4J 



\ ^1 ^2 ^3 ^4 / 



B 2J L 

Exactly as in |12j and |13j . the Central Limit Theorem follows from results in 
[42] . □ 

Lemma 4.8. Assume Condition]!] holds with < 4p + 2 — an < Q- Then, for 



S 2 log 2 S 



(b 2 iy 

Proof. The procedure is totally analogue to Lemma 19 in [13] , We obtain: 
Q M j{a) = Gj ,J^ )Gtj L (°)-(GLi(°))' 



(Gf ,J L (a)) 
1 K%(a)KV (a)-(KM(a)) 



2 
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/ Kf (,„) / / Kf (g) \ 2 Kf W) \ \ f Kf (op) 

I 2^3 JV J iff (a) I Z ^Af (a) y iff (a) I I JV J iff ( Q ) 

(V iv K v. (ao) f 



^ ; ' (^f(«)) 2 



Qj o Jl (a) can be rewritten as the sum of three terms: 

Qfo, J L (<*) = Qi (a) + Q2 (a) + Q 3 (a) , 



where: 



(a) 



)f° m (a) 



Kfjap) / Af l( a) \ / Af (a ) \ / Af M / Af l( a) 

^ JV ^«( a ) ^K-f(a) ^ l^j' V U«(«) y ^j iv J Af (a) ^ Af(a) 

2^3 JV J iff (a) J 

Q!T m (a) 



(E3 tf. 




Af(a) I Kf (q) ) I l^j JV J iff (a) I I l^j JV J I iff ( Q ) 



<f (a) y 
h (a) = ^ 



3 W 

den 



?2 



The next step consists in showing that: 

Q 2 (a) + Q 3 (a) = oj L (l) . 
Using Corollary [TJ Q2 um (a) can be written as: 

QT 



0™ m ( a ) 



^3 (") V 7 p.o( B ) VW*) 



■oj t (l) 



J 

while Q™ um (a) becomes 



3 (°0 
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so that: 



Q^ um (a) + Ql um (a) 

r\den . 



2 («) 

It remains to study Q^ n (a) ; by Propositions [1] and [51 we have: 



3 



lim 4^ ( S > « - «o) 



^2(1+^2-) ^2 

= 4 /p (B, a - ao )_ ( _^- y __ ? _>0. 

Finally, we prove that (c5_l) — > p TbT~[^ ■ Using again Proposition [T] and 
Corollary [TJ we write the numerator Q™ um (a) as: 



Kf (a) ') \f? ° Kf{a) V s ^/p,o(S) 



3 ■? 



Kf(a) V 6 IpAB) 



Let s = 2 (1 + S^a); by applying Corollary H we have: 

jJEoS*^* W = A^oo B« Zj - (S) 

= log 2 # ^ _ jT4 4/ (-B, a , a) . 



2 
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It remains to study Qf en (a) ; by using again (|A.6[) and Proposition [2] 



Jl— too 

B 2sJ L "«1 ^ ^sJx, 



= c B Jj, (B, a - a 
Hence 



73 s - 1 



B 2 ( 1 +^)log 2 i? 



2 



For the consistency of 3l, for o7l G [ao — Sl, ao + Sl], we have 

B 2 log 2 £ 



o,J L 'P , d o -, n2 



(52 _ i)- 



□ 



5. Narrow band estimates 

From Theorem 14.61 it is evident that, under Condition [31 the presence of the 
bias term does not allow asymptotic inference. As in [12] and [13] . we suggest 
a narrow-band strategy, developed only on the higher tail of the power spectrum, 
which allows us to avoid the problem due to the nuisance parameter. We start from 
the following 

Definition 5.1. The Narrow-Band Mexican Needlet Whittle estimator for the 
parameters fl — (a, G) is provided by 



Jl 



a,G 

j=Ji 



V B 2 N i / r-fl 



GKV in) J^ log \GK M 



V~, fc=l ^ 



or equivalently: 

(5-1) af lM = argmm< ;J >,G^ ;J >)), 

<;7» = ( log Gf 1 . jL (a)+ 1 f) iYilogi^f (a)) , 

where < Ji < Jl is chosen such that B Jl — B Jl — >• oo and 

(5.2) = B Ji (1 - fl (J,)) , Ji = Jl + l0g(1 - g R (JL)) . 

logis 

3 

We choose < g (Jl) < 1 s.t. Iimj^oo g (J L ) = and limj^oo J L g {Jl) 2 = . 

For notational simplicity B Jl is defined as an integer (if this isn't the case, 
modified arguments taking integer parts are completely trivial). For definiteness, 
we can take for instance g (Jl) — Jl 3 ■ 

Theorem 5.2. Letaj h -j x defined as in \5.1\) . Then under Condition\3\we have 
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2:S 



g (J L )i B Jl (a JvtJL -a Q )^N (o, ^y 1 ) 



, 2 B 2 ( 4 /logB-1 
<f> (£?) := log B 



where 

(B 2 -l) 2 \{B 2 -l) V logB 

Proof. The proof is very similar to the full band case, hence we provide here 
just the main differences. Consider: 



Qjf^t («) = ^K;J L («) ■ 

Let 

— , , -1 &K%(ao)*( /? 2 fe;p Gf ilJL (a o y 



S 



^ M M fcl V G ( Q o) *f («o) Go 
Simple calculations based on Proposition [5] lead to 



E**| = JW^ (B^-B^-V) +o(B 
\j=Ji 



(B 2 - l) n 
= B' Ul +0(B 2njL g(J L )) 



2nJ L \ 



We have: 



lim 



lim R,L ,. J P .o("o + l) 



J L g (J £ ) ip, (a ) 
E log^ ■ (S- - + ^)) j + o,, (1) 

Jl^oo i P) o(ao) -d + 1 

As in full band case, we collect out all the covariance terms defined in Lemma 
IB. 1 l and following Corollary [3J we have 

Var (Sj^j, K)) = - a ° (1 + T) N3 (.ZMr. (2) + 0j L (1)) . 
After some manipulations, we have: 



24 CLAUDIO DURASTANTI AND XIAOHONG LAN 

,2\ 2 ,\2 



(i-g(J L )f\ {B 2 -iy 2 2 

1 (l-ff(Ji)) (1 - log B (1 - g (J L ))) 



B 2 — 1 



B 2 J B 4 

2 



AZ jL . tJl {g{J L )) + o(g{J L f) 



^))g(J L) + 0(B^g(J L f) , 



where 

AZ Ji;Ji (//(J,)) = ( 1 + l"^ ) ~ ft ~ 2g (^)) ^ + i^gg < ^ ' 

4 VAl) f(Jl) 



(B 2 - 1) V iog- 8 , 

Thus we have 

Z Ji;Jl (2) = B 4J ^ (B) g (J L ) + O (s 4Ji 5 (Jl) 

Note that $ (B) > for B > 1. 
Hence we have 

far («o)) = ^o(l + ?)* ( B ) 9 (Jl) B^- , 

Consider now Qj l] j 1 (a), which we rewrite as 

Following a procedure similar to Lemma 14.81 we have 

Qnum (a) = c 2 B Gll p (B, a - a ) Zj u j l (s) , 
where s = 2 (l + 9 ^r iL )- Following ([531), we obtain 

Qnum (a) = c 2 B G 2 I (B, a - a ) $ (B, s) B 2s,,L g ( J L ) + O (b 2sJl 9 ( J l ) 2 ) , 
where 

$(5, S ) = log 2 5 ^(J.) , S logS-2 



(s s - 1) V fis - 1 lo & B 

Finally, we obtain 

Qden (a) = c|G2/ (B, a - ao) ( ^ B 8 - 7 

= c 2 B G 2 I p {B,a-a Q )B 2sjL +o(B 2sJl ] 

Hence 

Q^ ;Jl (a) = $ (-B, s) g (J L ) + O (b 2 ^ 9 ( J l ) 2 
and, for the consistency of a, we have 

Qjlj L (5)^ P <S>(B)g(J L ) ■ 
Thus _ 
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as claimed. Finally we can see that 



9 



(j L y b Jl e 



s 



Ji\Jl 



(«o) 



Q M 



Ji;Jl 



(a) 



O 



Jl- 



□ 



6. The Plug-in procedure 

In this Section, we will present a plug-in estimation procedure for the spec- 
tral parameter ao based on the interaction of the approach described here and the 
one based upon standard needlets introduced in [13!. As already mentioned in 
the Introduction, there already exist in literature Whittle-like estimators for spec- 
tral parameter based on spherical harmonics and standard needlets. The former, 
although characterized by a higher efficiency, can be affected by the presence of 
masked regions over the sphere, common set-up in Cosmological investigations, be- 
cause of the lack of localization in the spatial domain. The latter, as one here 
developed, is not altered by partially observed regions, paying the price of a lower 
precision. Therefore, our aim is to show that, if 4p— ao > , the spectral parameter 



M 

Jo, Jl 



estimator a 
aj L . First of all, observe that 



is more efficient with respect to the standard needlet estimator 



lim 

Jl— >« 



B 2 ' L Var 



M 

Jo, Jl 



a 



(B 2 



1 



Jl 



lim B 2jL Var (a,j L - a ) 



'B 2 



1 



B 4 log 2 

\3 



B 



B 4 log 2 B 



see again [13j. We can therefore observe that for 4p — ao > 0, 



2,2 

o-i < Po 



where a\ :— <t 2 (1 + t). Consider that, for any fixed p 
not depend on B 7 becomes, by the Stirling's formula, 

2 



4p > ao, ctq, which does 



V(2p-*20) 



We have that ctq is smaller than 1 for Ap > ao — 2, while easy calculations show 
that t < 1. On the other hand, as described in [13) . p\ — p\ (ao, B) is decreasing 
on B (see also Table 1): any attempt to reduce its value will produce an increase 
of the variance due to the term (B 2 — l) / B 4 log 2 B. 





Standard Needlet-pp 


Mexican Needlets - a\ 


B = $2 


B = V2 


B = 2 


p = 2 


p = 3 


p = 4 




= 2 


5.00 


2.24 


1.16 


0.62 


0.49 


0.42 




= 3 


5.04 


2.53 


1.34 


0.67 


0.51 


0.43 


a 


= 4 


5.10 


2.64 


1.57 


0.75 


0.55 


0.45 



Table 1: Comparison among different values of the variances pjj and Ctq. 
Hence, the plug-in procedure can be implemented in two steps: 

• First step: compute aj L in the standard needlet framework. 

• Second step: if p > a,/ L /4, compute a^ Jl by the mexican needlet ap- 
proach; otherwhise, accept aj L . 
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Appendix A. Auxiliary results: preliminaries 

The results collected in this section, provided by standard analytical calcula- 
tions, are here reported to explicit the structure and the behaviour of the function 
f p (-) defined in in (gig) . 



Lemma A.l. Let 



poo 

W 2a ,b,s = / t 2a exp (~bt 2 ) log s (t) dt 
Jo 



We have 



W 2a .b,o = 



&-(° + *). 



-r a 



w 2a 



b-( a+ i) r d 



6,1 



da 



logfe 



2 J ' 



r a 



&-(•+&) 



2a, b, 2 



_ 2 log b—+ log 2 b 
da z da 



T a 



Proof. Standard calculations lead to 

/>oo 

V^2a, ,o = / t 2a exp (~bt 2 ) dt 



b-( a +i) 
2 

&-(»+*) 



-r a , 
2 V 2 



(6i 2 )° 2 exp (-fa 2 ) (2fad*) 

r 



Similarly 

WW = 



/•OO 

/ t 2a exp (-bt 2 ) logtdt 
Jo 

l - £ t 2a exp (-bt 2 ) log (fa 2 ) dt- 1 ^ J™ t 2a exp (-fa 2 ) 



2 

b-(«+?) 
4 

6-( Q +i) 



c a 2 exp (—x) [log x — log b] dx 



— - log b 

da 



W a . 



b.2 



t 2a exp(-bt 2 ) {\ogtfdt 
b-( a +i) 



(bt 2 ) 2 exp (-fa 2 ) (log fa 2 ) - 2 log 6 log fa 2 + log 2 6 



Ibtdt 



_ 2 log 6— -flog 2 b 
da z da 



T a 



a 
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Lemma A. 2. Let f p (•) be defined as in \2.2]) . Then we have 

<-> E/;(i)-^r( O!)+ ^) +0 ( s * + ,) ; 



/>1 

Moreover, for Aj £ Z ; we have 



(A . 2) . *"»r m / a2 2+1 n o (s( „ +1)j) 



where 



„ dA, , „ R-Aj \ -((ai+a2)p+ : 4 i ) 



(Aj) , 

ai + a2 
Proof. Observe that 

/ 



l n 



e»p (-' 2 (tSj + -M) { ( w)T {{ ^i)T'- 



i>\ 



gjn f t l \2 / aiB 2A J +a2 \\ // j v 2\ («i+«2)P+7 

B 2o 2P Aj £ eX P J ^ ]^Aj JJ^bJJ J 

Rj'fn+l) / p2Aj \ (ai+a 2 )p+ I1 ± i .+00 

^ A . f — ^ ) / exp (-1) d* + o B^A 

2B 2a 2 pA 3 y aiB 2A 3 + a2 J J q V / 

B 2A,-(a lP +^±i) / + j X 
xrr (01 + a 2 ) P H +0 B (n+1) - 7 



1+1 ' p,ai,a 2 



(Aj)r(^+a 2 )p+^P) + o(i? (n+1)j 



2(ai + a 2 ) (Ql+a2)p+ 
Fixing Aj = 0, ai = a 2 = a/2, we obtain 

as claimed. □ 

We now investigate the behaviour of the function Kj 1 (a) and its derivatives. 
Proposition 1. Let 

V («) = 7^ WV+i-a.2. s = 77- / t 4 P +1 ~ a e- 2t2 (log i) s dt , s = 0, 1, 2 . 
Ob Ob jo 
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Then we have 

(A.3) Kf{a) = (I p , (a)+ Oj {l))B-^ 

2 -(2p-f+l) 



(A.4) = _(jlogB+^)+ (l))jrf (a) 



A -M (a): ^ Kf(Q) 



(A.5) = l^ log ^ + 2jlo g B^iM + iiM + <>(l) N )Kf (a) . 

Proof. These proofs follow the ones concerning the scalar needlet case (see 
We have indeed 

bw-l-y f-V p ex P (-2 (±y\ (m^ai 



«f («) 



i>i 



S- a %, (a) + Oj (S-^) , 



z>i 



C B B 2 3^\Bn * \ \Bi 



= -iff (a) logB J - -^B~ a] ^ P +i-a e -2t 2 bgMf + G ^ X ) 

«S«0 ^ ^fE(sr)*«*(->(sr) ,, )c«+')^V 

= Iff (a) (logB^) 2 + 2£- Q3 logS^ f ^- y t 4 P+i-« e - 2t2 \ogtdt + o(l) 



, 2 log 2 B + 2j log B^j + HI + W 
ip,o (a) ip,0 (a) 



□ 
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Corollary 1. From Proposition]^ we have that: 
(A.6) rWJ^) = !p ( B ' a ~ «o) i? (Q °-" b + o (flC«°-<*>Jj 

where 

I p (B,a-a ) := (2(2p+l))^ . 

Proof. The computation above shows that 

2 

Ob 

2-(2 P - f +l) 

.r(2p + i--) , 



and following (IA.3I) 

iff (a) , v ao-a T(2p+1 



Xf(oo) ' " r(2p+l-2a) 

as claimed. □ 



The next results follow strictly Proposition 27 in |13j . hence we will report the 
statements, while we will omit the proofs. 



Proposition 2. Let s > 0, B > 1. Then 
(A.7) £ B« ■ = (b sJi - - B'Vo-V) ; 



j=Jo 

(A.8) ±B»l\o g B> = ^— log B^J L -^- J )B^ 



((A-D-Bjb)^*- ') 



(A.9) J] B« (logfl') 



B s „ / / / 1 V S 



2 



Corollary 2. Lei 

vjosj* (*) = ( E /r; ] ( E /r/ o^/n 2 ) - ( E s^iogs^ 
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The we have 



v Jo ,j L 00 



B s log B 
B s - 1 

Moreover if Jq = — Jj, 



B s 



B sJ L _ B s(J -l)\ _ B s(J L + J -l) _ + jj2 



Vj L (s) : 
so that 
(A.10) 



B s log £? 
£ s - 1 



— — (V J - - B-(-A-l)V - J_^s(2Jt-l) (2J + x) 2 

5 s - 1 V J B s 



B 3s 



Urn B- 2 ' J *V> 1 (*)=log 2 J B- , 
Ji-H» (5 s - 1) 

Appendix B. Auxiliary Results: Covariance terms 



Lemma B.l. Letrs (Aj) 6e defined as in i2.8\) . Hence we have for 4p—a > 0, 
Jo < 0, 
(B.l) 



Jl 



Jl-3 



B 2 



So (Jl) := ]T fl 2 > £ r B (Aj) B a ^ = (1 + r ) B 2J - + o (B 2J <<) , 

j=Jo Aj—-J -j 

■Jl Jl —j 

Si ( Jl) ~ B ** lo § Bj E TB ( A ^) B<IQA ' ? 



(B.2) 



B 2 log S 
5 2 - 1 



(1 + T ) J L 



Aj=-J -j 
1 



5 2 - 1 



(1 + ri) £> + o (£? 



Jl Jl-3 

S 2 ( J L ) := J2 B2J 1o § 2 B ' Tb ( A ^) BQ0AJ 

3 = Jo Aj'=-,7o-j 



(B.3) 

S 2 log 2 B 
~ B 2 - 1 

where 



(5 2 - 1)- 



TO 



2(4p+l-a ) 
(_g(4p+2-a„) _ 1) ' 
( R4p+4-Q _ 1 

ri := 2 4p+1_ct ° ^ 



?2 •= 2 4p+1_Q ° 



(£4p+2-a _ iy 

W P (B) 



(B 4 p- 



ao+l _ 



and 



( B G B 4p- aa + 1) + S 4 B( 4 P- q o) ^3 B 4p-ao+l _ 6) + B 2 (a?( 4 P- q °) + l) + l) 



Moreover if we define 
Zj l 



B 2 + l 



S ( Jl) S 2 (J l ) - £ 2 ( J L ) , 
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we have 



where 



J L =0 (£2 - l) 4 



£ 2 

Proof. Let us call P = (4p + 1 — a ) and observe that: 

Jl-3 

y B a ° A >r B (Aj)-l 

Aj=J -j 

J L-j Jl-3 r>Aj 

Y 5 Aj >osh(A7io g mr p - 1 = 2 P V 

/ -1 -, Ji-J i 



P- 1 



- 2 P 



where we have considered the case Jq < . Hence we have, from Proposition [5] 

-1 -, j—Jo 

t ~ y B -(p+i)Aj 



1 fl _ R-(^+l)0'-./o) 

~ 5^+!) -IV 

while we have 

Jl-3 -. Jl-3 

y = ~ y b-( p -v a j 

= I (i_ B -(p-i)(Jl- 

B( p -V-l\ 

Consider now 



y B^r B (Aj) 

j=J Aj=-J - 



= E b2j + E B2j E 2 p i? A ^ p+i ) + y b 2 > y 2 f b- a « p -v 

j=Jo 3 = Jo Aj=J -j j=l Aj=l 

We have, given that P + 1 > 0, if J < 



y b 2 > y b a « p +v 

j = J Aj=J -j 
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-1 1 

n B (P+1) _ 1 \ L ) 

3= Jo 

- 1 ( ( R-2 _ d2(J -1)A _ - g(P+1) ' 7 " f d(1-P)(Jo-1) _ R-(l--P)^ 

- 1 - 1 V J B^-i-ll 

= (s 2J -) . 

On the other hand, 

Jl Jl -j 

£ B 2j B- Aj ' (p - 1J 

j=l Aj=l 



3 = Jo 



' - B2JL , B2JL ) + ° ( S2JL ) 



s( p -!) - 1 v^ 2 - 1 - r 

- ( (^■)(^- 1) ) E ' J ' t °^ 

Hence we have 

j=Jo Aj=—Jo—j 



B 2 -l \ (B p+1 — 1) 
= ^-jB^(l + ro) + o(B^) , 
Similar calculations lead to 

Jl Jl-3 

J^B^logB 1 J2 B Aj T B (Aj) 

j=Jo Aj=-J -j 

= jr B 2 > log Bi+2 P I £ B 2] log & ]T /^(P+i) + ^ bg B , 'jp b -a,-(p- 

j = Jo \j=Jo Aj=J -j j=l Aj=l 

where 

]Ti^ log B J ! ^b-Aj(p-i) 

J=l Aj=l 



7 = 1 V ' 



B( p -i)-l \(B*-l)(B p + 1 -l) JL (i?2_i) 2 ( B (P+D_i) 2 ' 

= ^iogs / i r!!!_!L) + (B2 j i} 
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while, if Jo < 

Y B 2 > log B j Yl B Aj(P+1) = o (B 2J ^) 

3 = Jo Aj=J -j 

We hence obtain 

jr,B 2 * log B* jr B A V B (Ai) 

j=Jo Aj=-J -j 



B 2 logBB 2J - ( jn — , 1 / 2 p (B 2 +( p + 1 )-l)\\ 



5 2 - 1 I v u/ B 2 - 1 
-B 2 log BB 2Jl ( , 1 



5 2 -l V 5 2 -l 
Furthermore, we have: 

Jl Jl-3 

Y B 2j log 2 B j Y b ^ t b (Aj) 

j = Jo Aj=-.T -j 

- ^ 5 2j log 2 fl' +2 P [ ^ B 2 ' log 2 B- 7 ^ B A * P+1 > + Y B 2j log 2 B- 7 ]T B~ A ^ p - 

3 = Jo \j=Ja Aj=Jo-j j = l Aj=l 

We observe that 

YB 2 nog 2 B> Y B- A ^ p -^ 

3 = 1 Aj=l 

1 / B 2 log 2 B„„.(f„ 1 \ 2 B 2 _\ b _ (p _ 1)Jl Jt 

1 ' ' i=i 



^ ^-^j + E^iog 2 ^ 

(B 2 -l) [ L [(B( p +V-l)J (B 2 -1) Jl [ 



(B( p +V - l) 2 



( B2 + \ 2 W ^ B) , )) + o(B^) 
\(B 2 - l) 2 - l) 3 J J V 7 



B 2 log 2 B / 9 ?i B 2 + 1 „\ / 9 7 \ 

= T^T) B ( To J! " 2 Ji + (F^ 2 j + ( ) 

where 

(B 6 /^- 1 (B( p -D + 1) + S 4 B( p - 1 ) (B 3 B P - 6) + B 2 (B^" 1 ) + l) + l) 
W P (iJ) = 

On the other hand we have 

Y B 2j log 2 B j Y S AJ(P+1) = o (B 2J -) , 

3 = Jo Aj=J -j 

so that 

Jl Jl-3 

Y B 2i !°g 2 B j Y B Aj T B (&j) 

3 = Jo Aj=-J -j 
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B 2 \og 2 B 97 (, _ x , 2(l+ri) B 2 + l , „A , . 

= Jb^T) B ( (1 + ro) Ji - W^Tf JL + W^f (1 + T2) j + ° ( ] 

Hence we have, from Corollary [2j that 

B 4 \oe 2 B A1 ( ( ' „ . t9 2(1 + ti) 5 2 + l , 

* = (^Tf 5 t ( v ( v (1 + To)J ^7^ Ji+ (^Tf (1 + T2) ' (1 + ro) 

i _ x2 

^(l+T0)--ga— j(H-Tl) 

5 6 log 2 i? 4/ / /(5 2 + l) „, 27x-7 2 \\ 

(5 2 -l) 4 

as claimed. □ 

Appendix C. Auxiliary Results: Cumulants 
Lemma C.l. Let Aj and Bj be defined as in and \4. 1 As — > oo 

-^cum | J2(A h + B n ), J2( A h + B^YSAh + B h)> + BjJ \ 

\ ll I2 ^3 ^4 ) 

f jj log 4 B 

Proof. It is readily checked (see also [12] ) that 

cum {6i,d u 61,61} = o(r 3 r 4ao ) . 

Let us compute: 



C 4 ,• „• = cum 



" • ' J " ' » G„#£ (00) ' W ia G Xjf (n,,) ' .V ;i G K% (a ) ' G (o ) 



h,h,hM \i=l Jl Js v ; / 

? (2 < + 1)4 (n % i^L) ) cum (°- o - o -°) + ( B_4i ) 



° fe fn^^^/p (J^)) s (2 - 4a °w (z 1 - 4 " )^ 
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Then we have 

'Gf o , jL (a ) G% tJL (a ) G^> ) G^>o) 



cum 



^0 <-*0 



o 



O 



( 77^7- E N nN j2 N ]3 N u Cf uj2j3j \ 
\ hhhu J 



As in [12] and |13j . the proof can be divided into 5 cases, corresponding respectively 
to 



I I . . . . I I 

-cum ■ 



B u L — 1 E A n » E A h > E A *> ' E A >* \ ' -^UI cum \ E B n > E S J2 . E fl i3 > E B i« 

jl h 33 34 I I 31 32 33 34 



B iJ 

and 



1 

cum 



E A ii ' E B n > E S i» ' E B « > , ^; c "™ ] £ . E 4fe ' E B h f E s « 

jl 32 33 34 ) y h 32 33 34 



1 

cum 



E A h ' E A ii ' E ^ 3 ' E 



Jl J2 33 34 

where we have used 14.101 14.111 We have for instance 



B i.j 



TJZ cum \Yl A ii > E A h ' E A h ' E A J' 4 f 

[ 31 32 33 34 ) 

\ 31,3233,34 »=1 / 

(^7T E l0 § 4 = °(^7T E V ^) = 



and 

TUZ cum I E S ii ' E S » ' E S ia . E B J4 



L\J2 B2J ^ 



31 32 33 34 

4 , ^ 

Gj L (a ) Gj L (a ) Gj L (a a ) Gj L (a ) 



cum 



B^\^ b Go ' Go ' Go ' Go 

3 

= O (log 4 B Jl B- 2Jl ) ; 
The proof for the remaining terms is entirely analogous, and hence omitted. □ 
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